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Abstract. We consider the Cauchy problem for a hyperbolic pscudodiffcr- 
cntial operator whose symbol is generalized, resembling a representative of a 
Colombeau generalized function. Such equations arise, for example, after a 
reduction-decoupling of second-order model systems of differential equations 
in seismology. We prove existence of a unique generalized solution under log- 
type growth conditions on the symbol, thereby extending known results for 
the case of differential operators f ll7H19l V 
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1. Introduction 

This paper establishes existence and uniqueness of a generalized solution to the 
scalar hyperbolic pseudodifferential Cauchy problem 

(1) d t u + A(t, x, D x )u = / when i <E (0,T), 

(2) u(0) = g. 

The data / and g are Colombeau generalized functions and A is a generalized 
pseudodifferential operator of order 1. Its symbol is represented by a family of 
smooth regularizations, which may (but need not) be convergent to a distributional 
symbol. Problem 11121) represents an extension of the scalar case of the partial 
differential equations considered by Lafon and Oberguggenberger in |17l 

One may think of problem (11121) as resulting from a system of second-order (par- 
tial differential) equations by reduction to first-order followed by a decoupling into 
scalar equations (cf. \22i Section IX. 1]). This is a standard technique in applica- 
tions, for example, in mathematical geophysics, where one decouples the modes of 
seismic propagation and subjects these to further refined analysis (cf. As they 

stand, these reduction-decoupling methods are rigorously applicable in the case of 
models with smooth coefficients or symbols, but cease to be well-defined under the 
realistic assumptions of only measurable (bounded) coefficients, which are to rep- 
resent the elastic or acoustic properties of the earth's subsurface. Moreover, the 
initial value and the right-hand side are distributions corresponding to the origi- 
nal seismic source and force terms, which are, by nature, strongly singular, e.g., 
delta-like. If the original model coefficients are replaced by regularizations, then we 
may carry out all transformations within algebras of generalized functions from the 
outset and arrive at (11121) in a well-defined way. The purpose of the current paper 
is to investigate the general feasibility of rigorously solving the resulting decoupled, 
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so-called one-way wave equation, by generalized functions. Future work will be 
devoted to the regularity analysis of the solutions and their asymptotic relations 
with distributions. 

A word on conventions and notations concerning the Fourier transform: if u is a 
temperate distribution on R" we denote its Fourier transform by u or Ju; occasion- 
ally, when several variables and parameters are involved, we write expressions of 
the form J x ^^{u{y, x)) to indicate that the transform acts on the partial function 
(or distribution) u{y,.) and yields a function (or distribution) in (?/,£); the inte- 
gral formulas for the transforms follow the convention Ju{— x)/(2ir) n — 3 r ~ 1 u(x) — 
J exp(ia;£) u(f) <f£, where cf£ = d£/(27r) n . 

Subsections 1.1-3 serve to review Colombeau theory, fix our notations for general- 
ized symbols, and also recall the corresponding result on the Cauchy problem for 
hyperbolic differential operators with generalized coefficients. Section 2 establishes 
precise energy estimates, which are at the heart of the existence and uniqueness 
proof for the Cauchy problem presented in Section 3. Finally, under additional as- 
sumptions on the symbol and data regularity, we are able to draw some conclusions 
about the solution regularity which are revealed by the technique of the existence 
proof itself. 

1.1. Colombeau algebras of generalized functions. We will set up and solve 
the problem in the framework of algebras of generalized functions introduced by 
Colombeau in E] . More specifically, we will work in a variant which is based 
on L 2 -norm estimates as introduced in pQ. We will recall the definition and basic 
properties below. As general references and for discussions of the overall properties 
of Colombeau algebras we refer to the literature (e.g. |31l51l20p. 

We consider the space-time domain At '■= R™ x (0,T). The basic objects defining 
our generalized functions are regularizing families (u e ) eg (o,i] of smooth functions 
u e E H°°(Xt) for < e < 1, where H°° denotes the intersection over all Sobolev 
spaces. To simplify the notation, we shall write (u £ ) £ in place of (M £ ) e g(o,ii through- 
out. We single out the following subalgebras: 

Moderate families, denoted by £m,t 2 (At), are defined by the property: 

Va € Nq 3p > : \\d a u s \\ L2 = 0{e- p ) as e -> . 

Null families, denoted by N^At), are the families in £m,l 2 (^t) having the fol- 
lowing additional property: 

Vq > : \\u £ \\ L2 = 0{e q ) ase^O. 

Hence moderate families satisfy L 2 -estimates with at most polynomial divergence 
as e — > 0, together with all derivatives, while null families vanish faster than any 
power of £ in the L 2 -norms. For the latter, one can show that, cquivalcntly, all 
derivatives satisfy estimates of the same kind (cf. [Hj)- The null families form a 
differential ideal in the collection of moderate families. The Colombeau algebra 
Sl 2 (Xt) is the factor algebra 

S L 2(X T ) = £. M J y 2(X T )/>{ L 2(X T ) . 

(The notation in pQ is 82,2, and correspondingly for moderate and negligible nets, 
where the variability of L 9 -norms in the definitions was essential.) The algebra 
Si2(]R n ) is defined in exactly the same way and its elements can be considered as 
elements of Sl 2 (At)- On the other hand, as explained in Remark 2.2(i) and 
Definition 2.8], the restriction of a generalized function from Sl 2 (At) to t = is 
well-defined: for any representative {u £ ) e in E ML 2(Xt) we have u £ € C 00 (R™ x 
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[0, T]) (i.e., smoothness up to the boundary of the time interval) and that (u e (., 0)) e 
belongs to £M,L 2 (K n )' We use the bracket notation [ . ] to denote the equivalence 
class in Sl 2 ■ 

Distributions in H-°°(W l ) = \J seR H s (W l ) are embedded in 9 L 2{R n ) by convolu- 
tion: l(w) = [(w * (p e )) £ ], where 

(3) Pe {x)=e- n p(x/e) 

is obtained by scaling the fixed mollifier p, i.e., a test function p G J^(R n ) of integral 
one with all moments (of order 1 and higher) vanishing. This embedding renders 
#°°(R n ) a faithful subalgebra (cf. Theorem 2.7(h)]). In fact, given / G H°°(W l ), 
one can define the corresponding element of Sl2(R' 1 ) by [(/)J (with representative 
independent of e). In the same way we may consider H°°(Xt) a faithful subalgebra 
of Sl 2 {Xt)- 

Some Colombeau generalized functions behave macroscopically like a distribution. 
We say that u = [(u E ) e ] G 9l 2 is associated with the distribution w G CD', denoted 
by u ~ w, if u e — > w in D' as e — > 0. 

Intrinsic regularity theory for Colombeau generalized functions has been a subject 
of active research. Its foundation is [2111 Section 25] with the definition of the 
subalgebra S°° of S, which plays the same role for S as C°° does within D' . The 
basic idea is to couple the generalized regularity notion to uniform e-growth in all 
derivatives and it leads to the important compatibility relation 

3°° nay = C°°. 

Similarly, we define here the subalgebra S^°2 of regular elements of Sl 2 by the 
following condition: u = [(u e ) e ] G 5l 2 belongs to 3^2 if and only if 

(4) Bp > OVa G N' l : \\d a u e \\ L2 = 0(e- p ) as e -> . 

Observe that p can be chosen uniformly over all a. In particular, if u £ — v * p e with 
v G H°°, then p = is possible when we let fall all derivatives on the factor v. 

Concerning sources for recent and related research in Colombeau theory, with a 
diversity of directions, including such topics as pseudodifferential operators with 
generalized symbols, regularity theory, and microlocal analysis of nonlinear singu- 
larity propagation we refer to |71 ITT1 IH1 IH1 IHl ITS] . 

1.2. Generalized pseudodifferential operators. For comprehensive theories of 
approaches to pscudodiffcrcntial operators with Colombeau generalized functions 
as symbols we may refer to the recent literature on the subject [7||H1EI- However, 
the purpose of the present paper is to present a short and self-contained discussion 
of the solution to the hyperbolic pseudodifferential Cauchy problem. Therefore 
we do not need to call on the full theory of generalized symbol classes, mapping 
properties, and symbol calculus, as it has been extended systematically and with 
strong results in Nevertheless, this background will be substantial in fur- 

ther development, refinements, and applications of the current work, in particular, 
concerning regularity theory and microlocal analysis. 

We will use families of smooth symbols satisfying uniform estimates with respect 
to the x (and t) variable as described in |1(JM16) . To fix notation, let us review the 
definition. A complex valued function a G C°°(1R™ x R") belongs to the symbol 
class S m of order m G R if for all (a, (3) G Nq™ 

(5) C»:= sup (l + |e|)- m+ l Q '|^a(x,OI <oo. 

(a:,f)eK»+» 
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S m is a Frechet space when equipped with the semi-norms 

(6) g£i(a):= max c^(a), 

|a|<fe,|p|<i 

a notation we will make use of freely in several estimates in the sequel. (Ob- 
serve that compared to the semi-norms and notation used in we have = 
max{q™ r (a) : k + r < I}.) In fact, we will use symbols which depend smoothly on 
time, considered as a parameter. More precisely, we consider the space of symbols 
a(t,x,£_) where a € C°°([0, T], S m ) (i.e., each ^-derivative on (0, T) is continuous 
into S m up to the boundary t — and t — T) with the semi-norms 

(7) Q™ fc >):= max sup •=•))• 

0<i<3 ts [o,T] 

By a generalized symbol we mean a family (a £ ) ee [ 0j i) of smooth symbols in S m (the 
same m for all e) which satisfy moderate semi- norm estimates, i.e., for all k and I 
in N there is N G N such that 

(8) 0( £ -^) ( e -0). 

Generalized symbols with parameter t G [0,T] are given by families (t,x,£) i— > 
a e (i,:r,£) (e G (0,1]) such that a e G C°°([0, T], S m ) with moderate semi-norm 
estimates: for all j, k, I in No there is N G No such that 

(9) Q% ktl (a s ) = 0(e- N ) (e - 0). 

Obviously, no major changes would be required to incorporate more general types 
of symbols, especially the Hormander's classes S™ s would mainly require changes 
in notation (at least when < p < 6 < 1). 

Let (a e ) e be a generalized symbol with parameter t G [0,T]. We define the corre- 
sponding linear operator 

A : 9l 2 (Xt) — > Sl 2 (Xt) 

in the following way. On the representative level, A acts as the diagonal operator 

(u E ) e i ^ (a e (i, a;,L> x )-u e ) e V(u £ ) £ G Sjw.i^Xr). 

Here, a e (t, x, D x ) acts as an operator in the a; variable with parameter t. The 
moderateness of (a £ (t, x, D x )u e ) e follows from JHJ) and the fact that operator norms 
of d\dP o a e (x,t,D x ) on Sobolev spaces are bounded (linearly) by finitely many 
semi-norms of the symbol (cf. [TBI Ch. 3, Theorem 2.7]). In the same way, it 
follows that null families are mapped into null families, so that A is well-defined 
on equivalence classes. We call A the generalized pseudodifferential operator with 
generalized symbol (a E ) E . 

1.3. Review of hyperbolic partial differential equations with generalized 
coefficients. We briefly review the situation for symmetric hyperbolic systems of 
partial differential equations in Colombeau algebras. The heart of this theory was 
developed in |17U19| . from where we recall the main result on the Cauchy problem. 

The theory is placed in S instead of Sl 2 , i-e., the data /, g as well as all coefficients 
satisfy asymptotic local L°°-estimates of the kind described in the introduction. In 
view of our intended generalization of the scalar case to pseudodifferential operators, 
let us simply focus on this situation in the Cauchy problem (11121 . We have / G 
S(R n+1 ) and g G S(R") and the spatial operator A is a differential operator of the 
form 

n 

A = ^a ] (x,t)d X3 +b(x,t) 
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where the coefficients dj, b are in S(IR n+1 ), a,j real. Note that a generalized symbol 
for A is given by 



where dj lS , b £ are any representatives of aj, b; a^ e may taken to be real- valued. 

Sufficient conditions for existence and uniqueness of a solution u £ S(M n+1 ) to illlL'l) 
are as follows: 

(i) aj, b are equal to a (classical) constant for large |jc| (any kind of uniform 
boundedness in x and e for large \x\ would do; it ensures uniqueness and 
enables one to use partition of unity arguments in the proof) 

(ii) b as well as D^cij are of log-type, i.e., the asymptotic norm estimates (of 
order 0) have bounds 0(log(l/e)) (this ensures existence by guaranteeing 
moderateness from energy estimates). 

Counter examples show that none of the two conditions can be dropped without 
losing existence or uniqueness in general. 

Remark 1.1. It turns out that the non- uniqueness effect as constructed in |191 
Example 1.4] disappears in Sl 2 - (I n mentioned example, the constructed initial 
value v{x,0) = [{x{x + 1/e))], X e 3)(E) with x(0) = 1, is negligible in S but 
gives \\v e (., 0)\\ L 2 = ||xllz,2 > 0, hence is nonzero on Sl 2 -) As a matter of fact, 
the L 2 -energy estimates, to be discussed in the following section, directly yield 
uniqueness; this holds even with coefficients that allow for logarithmic growth as 
e — > throughout the entire domain. 

The non-locality of pseudodiffcrential operators seems to prohibit an adaption of 
the proof technique of |17| . where one is able to pass from L 2 -energy estimates 
to local L°°-estimates. On the other hand, when working in 5l 2 j there is also 
the structural advantage of having good mapping properties of pseudodifferential 
operators with uniform symbol estimates on Sobolev spaces. 



Our proof of unique solvability of the Cauchy problem will be based on energy 
estimates, with precise growth estimates of all appearing constants depending on 
the regularization parameter e as e —> 0. This in turn is solely encoded into the 
generalized symbol in form of the semi-norm estimates of the regularizing (resp. 
defining) family of symbols. Therefore, and also to make the structure more trans- 
parent, we will first state the preparatory estimates for smooths symbols in terms 
of explicit dependencies on symbol semi-norms and insert the e-asymptotics only 
later on. 

In order to maintain close resemblance in notation with the cases of differential 
operators or decoupled systems, we shall write the symbol of A in the form i a(t, x, £) 
with a £ C°°([0, T], S 1 ); in other words, we review energy estimates for the operator 



n 



(10) 




2. Preparatory energy estimates 




P := d t + ia(t,x,D x ) 



with eti real-valued, clq of order 0, 
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or equivalently, that 



{I2J' 



a{t,x,D x ) - a(t,x,D x ) 



is of order 0. 



Besides stating the general case in the following proposition we also give details on 
two special instances. These are of interest in applications and allow for certain 
improvements concerning the regularity assumptions in terms of symbol derivatives, 
which are required in the constants of the basic energy estimate. 

Proposition 2.1. Assume that P is the operator given in 1)11(1 and such that 
G3 holds. Let u e C([0, T],H 1 (W n )) n C 1 ([0,T],£ 2 (K™)) and define / := Pu £ 
C([0,T],L 2 (R™)). Then we have the energy estimate 



where the constant C > as well as k' n , l' n , k n , l n are independent of u and can be 
chosen according to certain assumptions on the symbol a as follows: 

(a) General case: We have k' n = l' n = [n/2\ + 1, k n = i(\n/2\ +1), l n = 
2{n + 2) and C depends only on the dimension n. 

(b) Constant for large |x|: If there is ro > such that 

(14) a(t, x, £) = h(t,£) whenever \x\ > ro, 

where h is a symbol of order 1 (with parameter t and no x variable) , then 
C depends only on n, ro and the semi-norm orders are at most k' n = 0, 
l' n =n + l,k n = 1, l n = n + 2. 

(c) Real symbol: If in addition (to any of the assumptions above) the sym- 
bol a is real- valued, so that ao is real as well in 1)12(1. then the term 
Qo,k' .1' ( a o) can be dropped in (fbl|l . 

Proof. Using the standard decomposition of the operator a\(t, x, D x ) into self- and 
skew-adjoint part, a\ = (a\ + a\ )/2 + (ai — ai )/2, we obtain 



= 2Re{f(t)\u(t)) +21m{a 1 (t,x,D x )u(t)\u(t)) +21m(a (t,x,D x )u(t)\u(t)) 
< \m\\l> + \Mt)\\v + \\ai{t,x,D x )-a x {t,x,D x f\\ \\u{t)\\ 2 L2 

+ 2\\a Q {x,D)\\\\u{t)\\ 2 L2 . 



The operator norms (taken with respect to L 2 ) are finite by 112(1 and 1(12(1 ' (cf. (101 
Theorem 18.1.11] or |lfil Ch.2, Theorem 4.1]) and we will derive explicit estimates 
for these. 

For the proof of case (c), observe that we have Im(ao(i, x, D x )u(t)\u(t)) = if ao 
is real, so that the last term on te right-hand side of ((15(1 can be dropped from all 
further considerations. 

Case (a): We use a representation of the symbol of bi(t,x,D x ) :— ai(t,x,D x ) — 
ai(t, x, D x ) with integral remainder terms as it is developed in |161 Ch.2, 1-3] or 
[SJ Ch. 1,5-6]. According to this (or as sketched in the Appendix below), the zero 



(13) 




(15) ±\\u(t)\\ 2 L2 =2Rc(d t u(t)\u(t)) 
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order symbol b\ (t, x, £) is given by 



(16) bi(t,x,£) := ai(t,x,£) - ai(t,x,£) - I rj t e(t,x,(,)d6 

3=1 ■ 



n ~i 
3=1 J ° 



since di(i, a:, £) is real-valued, where 



(17) r i , fl (t,a: ) 0= / / e"*""" ^fl Xj oi(t,a; + y,£ + 9rj) dyctr, 
in the sense of oscillatory integrals. 

As a close inspection of the proof of ^3 Ch.2, Lemma 2.4] shows (which we detail 
in the Appendix), we have the following estimate for all d e No, (a, (3) S N§" 

(18) l^^r^C^^O] < Cd,^ (1 + |C!)- |o:| QWa+lcl^+a+ic.i+^iCai), 

which is uniform with respect to 9 € [0, 1]. Combined with formula l|16|) for b\ this 
yields 

Qt),fc,z(frl + 2 °o) ^ 2< 9o,fc,i( a o) + Cfc,i ( ?0,n+2+fc,n+2+fe+i( a l)- 

By the theorem of Calderon-Vaillancourt (or one of its variants, cf. Ch.I, Theoreme 
3], ^H], |H Ch.3, Corollary 1.3]), we have the general L 2 -operator norm estimate 

||6(t,X,Da;)|| < C I 'i ( 9o,Ln/2J+l,Ln/2J+l(^) 

whenver b € C°°([0, T],S°). Therefore we conclude that 

||6i(t, x, D x )\\ + 2\\a (x, D)\\ < C n (<2o,L™/2j+u«/2j+i( a o) + Qo,w( a i), 

for any k > n + [n/2\ + 3 and I > n + 2[n/2\ +4. This completes the proof of the 
general case. 

Case (b): Let x G £>(M") such that x( x ) = 1 for M < r o and < \ < 1- Then the 
term a(£, cc, £) — a(i, x, £) occurring in 116|l can be written in the form 



X(z)(a (i,x,0 - ao(t,a;,0) + (1 - - &(*,0) 

:=6o(*,a;,0 + (l-x(a:))/io(t,0- 
The second part in this decomposition is the operator symbol of a convolution with 
bounded symbol (since := h — h is of order 0), composed with multiplication 
by 1 — X from the left. Hence the L 2 operator norm corresponding to this second 
summand has the following upper bound 

\\(l- X (x))h (t,D x )\\ < II 1 - XIL- INLoo < 2||ao|| £ c 

Note that bo is a symbol of order zero with support contained in \x\ < ro. We 
will estimate the operator norm of Bq := bo(t,x, D x ) (on L 2 (M"), uniformly with 
respect to t G [0, T]) via the Schwartz kernel Kq of the "Fourier transformed" 
operator Bq := JoB o J^ 1 and using the fact that 

||B || = (2^r/ 2 ||B ||. 

As a distribution in C°°([0, T], y'(M. 2n )), the kernel is computed from the symbol 
by the formula 

(19) K (t, & r,) = (2tt)-" J(b (t, ., tj)) (£ - r,). 

Since x >— > bo(t, ., rf) has compact support it follows that -Ko is smooth on [0, T]xl 2n ; 
in fact, we will sec that it is an integrable kernel and hence we may apply a classical 
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lemma of Schur (cf. |101 Lemma 18.1.12]). Before doing so, we will first show that 
the remainder terms in (jlfij) are of a similar form. 

Consider formula 1)17(1 and introduce the short-hand notation bj := d^.d x M. Then 
z i — * bj{t, z, C) has compact support in \z\ < ro and we may write 



: ,e(t, x, = ?~l x (?(bj(t, + 6r,)){r,) 



Now let Rj_e '■= fj,e{t, x, D x ) and define, exactly as above, the corresponding oper- 
ator Rj£ with intertwining Fourier transforms; denote by Kj t g its Schwartz kernel. 
The above representation for the symbol r^g in terms of bj and direct computation 
yields the formula 

(20) K^ e (t, £, v) = (27T)-" ?(bj(t, ., 77 + 0(£ - »?))) (£ - v)- 

Equations 1(19(1 and 1(20(1 have the following structure in common: we have a symbol 
d G C°°([0, T], S°) which vanishes when |x| > ro and a smooth kernel K defined by 

K(t, £, rj) ■= (2tt)- b J(d(i, ., /(£, ??))) K ~ »?), 

where / : R 2 ™ — ► K™ is a linear map. In order to apply Schur's lemma we estimate 
the partial L -norms of the kernel and obtain 

(27T)" / \K(t,£,ri)\d£= [ |?(d(i, .,/(£, 77))) (£-r« 



and similarly 



= / \7{d(t,.J^ + rj, V )))(0\d^< I sup|y(d(t,.,C))(0|de. 



(2tt)" / \K(t^,r,)\dr,< \ sup \?(d(t, ., C))(r?)| dry. 



Assertion: There exists a constant c(n, To), depending only on n and ro, such that 
(21) f S up\?(d(t,.,0)(»)\d(i<c(n,r )Q o in+lfi (d). 

R* 

Noting that / \d£d(t, x,Q\dx < c n r^ \\d^d(t, ., C)Hl<x>> the proof is exactly as in [THl 
Theorem 18.1.11']. 

In summary, applying 1(21(1 and the general integral kernel estimates above to the 
kernels given by ((19(1 and ((20(1 (note that bj involves first-order derivatives of a in 
x and £ already) we have proved the claims of case (b) in the proposition. □ 



3. COLOMBEAU SOLUTIONS 

We return to the scalar pseudodifferential Cauchy problem 

(22) d t u + Au = f inX T , 

(23) u(0) = g, 

where X T = R" x (0,T) and with data / G 9l'(X t ) and 5 G S^(K"). A is a 
generalized pseudodifferential operator of order 1. More precisely, we assume that 

A: Q L 2(X T ) -*• S L 2(X T ) is given by (u e ) s i-> (i a £ (t,x, D x )u e ) ei 
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where (a £ (t, x, £)) e is a generalized symbol of order 1 with parameter t E [0, T]. In 
addition, we impose the hyperbolicity assumption 

(24) Ve E (0, 1] : a £ = eti j£ + ao i£ ai j£ real-valued, ao l£ of order 0. 

The semi- norms in the basic energy estimate (|13[) now depend on e E (0, 1], and, 
upon applying Gronwall's inequality, will appear as exponents in the L 2 -norm es- 
timates of a prospective generalized solution; this suggests to assume logarithmic 
bounds on the symbols. We say that a generalized symbol (fe £ ) £ of order m (with 
parameter t E [0,T]) is of log-type up to order (k,l) if 

(25) QfcM = 0(log(l/e)) (e-0). 

Theorem 3.1. Let A be a generalized first-order pseudodifferential operator, de- 
fined by the generalized symbol (za £ ) £6 ( 0i i] with parameter t E [0, T], and satisfying 
the hyperbolicity assumption Assume that (ai, £ ) £ is of log-type up to order 

(fc n , l n + 1) and that (ao, £ ) £ is of log-type up to order (k' n , l' n ). 

Then for any given / E Sl 2 (^t), 5 G Sl 2 QR") the Cauchy problem I|22l23jl has a 
unique solution u E S L 2(X T ) if fc„ = 3([_n/2j+l), l n = 2(n+2), = /; = Ln/2j+l 

Furthermore, we have variants of the log-type requirements in the following two 
cases: 

(i) If there is ro > and an x-independent generalized symbol (h e (t, £)) £ such 
that 

(26) a E (t,x,(i)=he(t,g) when |x| > r , 

then we may put k n = 1, l n = n + 2, fc^ = 0, l' n = n + 1. 

(ii) If a s is real- valued for every e E (0, 1] then no log-type assumption on ao, £ 
is required. 

Proof. Let (g £ ) £ E 3, (/ £ ) £ E / be representatives. At fixed, but arbitrary, s E (0, 1] 
we consider the smooth Cauchy problem 

(27) d t u £ + ia £ (t,x, D x )u £ = f s in X T , 

(28) u £ (0)=9e- 

It has a unique solution u £ E C°°([0, T], H°°(W 1 )), thus constituting a solution 
candidate (u e ) e (cf. QII Ch.7, Theorem 3.2] or Ch.6, Theorem 2.1] with ad- 
ditional i-regularity following directly from the equation). We have to show that 

(U e )e € £m,L*(Xt)- 

Denote by C £ := C(l + Qq k , v (ao, £ ) + k n ;„( a e)) tne constant occurring in the 
energy estimate (|13f) applied to u £ . Gronwall's lemma implies 

(29) \\u £ (t)\\l 2 < (\\g £ f L2 + 1^ \\f £ (r)\\l 2 dr) exp(C £ T). 

By hypothesis we have C £ = 0(log(l/e)) as e —* 0. Thus we obtain uniqueness 
immediately from (|29(l - once moderateness is established - because null family 
estimates for f £ , g £ then imply such for u £ as well. 

For the proof of existence, we first observe that the basic estimate for ||w £ || L 2( Xt ) < 
Tsup tg [ ,T] II^eWIIls = 0{e~ N ) follows at once from (^Hjl by the moderateness of 
the data. It remains to prove moderateness estimates for the higher order deriva- 
tives of u £ . 
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x- derivatives: Let 0^a£ Ng and apply d x to equation (|27|1 . It follows by induc- 
tion and simple commutator relations of a £ (t, x, D x ) with d x . that this produces 
an equation of the following structure. Denote by e 3 - = (Sj,k)k=i the j th standard 
basis vector in K™ then 

(30) d t d^u e + ia e (t,x,D x )d^u s + i (d Xj a lte )(t,x,D x )d^- e iu e = F eta , 

l<j<n 

where F e>a equals the sum of d x f s plus, if |a| > 2, a linear combination of terms 
of the form 

(31) (<9f a e )(i, x, D x )d^- u £ with < a and 2 < |/3|, 
and 

(32) (a xj ao je )(*, z, D x )d x ~ e] u £ where a,- ^ 0. 

Assume that moderateness of ||93'M e || i 2 has been established already when |— y | < |a|. 
Since <9f? a £ is of order 1 we have 

\\(dSa e )(t,x,D x )d^u e (t)\\ L2 < CiQS iTn>ro (a£a E ) 

< Cl Qo,m,m( 9 x a e) . max ||<9> £ (t) || L 2 , 

|7|<|q| 

where Ci, C(, and to depend only on the dimension n f |161 Ch.3, Theorem 2.7]). 
Similarly, since d Xj ao^ £ is of order we also have 

\\(d Xj a , e )(t,x,D x )dg- e 'u e (t)\\v < C2Q° , m/ , m ,(d Xj a 0tE ) \\dr^u £ (t)\\ L2 

< C^Qo.m'.m'idx.ao^) max ||<9> £ (i)|| L2 , 

where C'2, C 2 , depend only on the dimension and m! = [n/2j + 1. Hence, by 
the induction hypothesis, we have H-F^aWH/^ = 0(e~ N ) as e — > uniformly in 
t e [0, T] for some N. 

We return to equation l|3()|) . consider it as an equation for v e := d x u £ , and supply 
the initial value v £ (0) = d x u £ (0) = d x g £ . Applying the basic technique from the 
beginning of the proof of the energy estimate to equation (|3U|) we obtain 

j t \\v e (t)\\^=2Re(d t v E (t)\ve(t)) 

< H^aWII^ + K(*)ll£* + \\a £ (t,x,D x ) - a s (t,x,D x )*\\\\v s (t)\\ 2 L2 
+ E - (^«l t e)(t,^, ,D a; )*|[([[^- e ^ E (t)||i 2 + n^(t)iii 2 )- 

l<j<n 

If we define G(d x a 1>e )[t) := Y!j=i \\{9 Xj ai t£ ){t,x,D x ) - (d Xj ax, e ){t, x,D x )*\\ then 
we get 

Y f Pe(t)\\l < \\F s , a (t)f L2 +G{d x a 1<£ )(t) max \\d2u e (t)\\l* 

at It 1 < 1 a 1 

+ (l + ||<z £ (t, x, D x ) - a £ (t, x, D x f\\ + G(d x a 1<e )(t)) \\v s (t)f L2 . 
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The term H s>a (t) := \\F e , a (t)\\ L3 + G(d x ai, s )(t) ■ max| 7 |<| a | \\d2u e {t)\\ L2 is of mod- 
erate growth and, by the proof of Proposition l2.il we have 

1 + \\a e (t, x, D x ) - a e (t, x, D x f\\ + G(d x a h£ ){t) 

n 

< C (1 + Q° Wn (a , s ) + Ql, kn , ln (a s ) + Y,Q 1 oM n (d Xj a h£ )) := C e , 

j=i 

which is a log-type constant by the hypotheses of the theorem. Note that the 
specifications of k n , l n , k' n , l' n for the general case match those of Proposition l2.il 
case (a), whereas the hypotheses in (i), (ii) match cases (b), (c) there. Thus, we 
prove all assertions of the theorem simultaneously when the notation is understood 
in this way. Finally, integration with respect to t and Gronwall's lemma yield the 
estimate 

\\d>e(t)\\ 2 L2 < (\\d^g\\ 2 L2 +J\ e . a (r)dr) exp(^T) = 0(e- M ) 

for some M and e sufficiently small. Hence \\d x u £ \\ L2 ^ XT ^ satisfies a similar estimate. 
In particular, we have the same bounds on the spatial Sobolev norms ||u e (i)|| fffc for 
k arbitrary and uniformly in t £ [0, T]. 

t- and mixed derivatives: Equations Q27JI and (|30[1 directly imply estimates of the 
form ||0 t <9"it e (i)|| i2 = 0(e~ N ) for any a € Nq (and uniformly in t). To proceed to 
higher order i-derivatives, we simply differentiate equations l|27|l . resp. (|3U[) . with 
respect to t. The Sobolev mapping properties of the operators (dtd@)a e (t, x, D x ) 
and moderateness assumptions on the symbols then yield the desired estimates for 
\\d l t d x u e {t)\\ L 2 successively for I = 0, 1, 2, . . . and a arbitrary. □ 

Remark 3.2. (i) The key assumptions in Theorem l3.1l are the log-type estimates 
on the symbol. We know already from the differential operator case that they 
cannot be dropped completely. However, these are sufficient conditions and merely 
reflect the various operator norm bounds available for zero order symbols (as used 
in proving the energy estimates). Thus they cannot be expected to be sharp. In 
fact, the value of the theorem lies in a general feasibility proof and any special 
structure inherent in a concrete symbol under consideration in applications might 
allow for improvement. 

(ii) In order to meet the log-type conditions of the above theorem in a specific 
symbol regularization one may call on a re-scaled mollification as described in |19|. 
To illustrate this procedure, let us assume that the non-smooth symbol of order m 
is given as the measurable bounded function a(x, £) such that for almost all x the 
partial function £ i— > a(x, £) is smooth and satisfies for all a S Ng an estimate 

\\d^a(.,0\\ LOO <C Q (l + |£|r-K 

Let p be a mollifier and let < uj £ < (log(l/e)) 1 /' c for some k £ N, uj e — » oo as 
e — > 0. Let p E (y) '■= Lo £ p{u) e y) and define the regularized symbol by a e (x,^) := 
{p e * a 6 (., £,))(x) (convolution with respect to the x- variable only). Then it is easy 
to check that a e S S m and of log-type up to order (oo, k). 

As in [Ej j essentially by inspection of the above existence proof, we establish com- 
patibility with distributional or smooth solutions, that is macroscopic regularity in 
a certain sense, when the symbol is smooth. 

Corollary 3.3. In Theorem 13.11 assume that A is given by a smooth symbol 
a G C oo ([0,T],S' 1 ), i.e., a £ = a for all e £ (0,1]. 
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(i) If / e H°°(X T ), g £ H°°{R n ), then the generalized solution u £ Sl 2 (X t ) 
is equal to the classical smooth solution. 

(ii) Let / £ C([0,T],H s (R n )) and g £ H s (R n ) for some s £ M, and v 
be the unique distributional solution to (111211 in C([0, T], H s ). Define 
generalized data for problem I|22I23[) by / := [(/ E ) e ] £ 5l 2 (Xt) (resp. 
g := [(g e )e\ £ Sl»(R")), where f £ £ H°°(X T ) (resp. g £ £ H°°(R n )) are 
moderate regularizations such that f £ f in C([0, T], i/ s (R")) (resp. 
g £ — > <? in H S (W 1 )) as e — > 0. If u is the corresponding generalized solution 
in Sl 2 (Xt) then it is associated with the distributional solution v. 

Proof. Part (i): Since we may choose the constant nets (f) £ , (g) £ as representatives 
of the classes of / and g in Sl 2 , and a £ — a by assumption, we obtain the classical 
smooth solution to equation H27I28[) as a representative of the unique Colombeau 
solution. 

Part (ii): The unique solution v £ C([0, T],H S ) to 1112(1 depends continuously on 
the data / and g by the closed graph theorem. Hence the solution representative 
u £ , defined as the solution to (|27I28I) . converges to v in C([0, T], H s ) as e — > 0. □ 

Finally, we prove that the intrinsic regularity property for the generalized solution 
holds if the data are in S^ 2 and the generalized symbol is only mildly generalized, 
namely satisfies additional slow scale conditions. This notion was introduced and 
investigated in some detail in and found to be crucial for regularity theory of 
partial differential equations. Recall that a net (r e ) e of complex numbers is said to 
be of slow scale if it satisfies 

Vp>0: \r e \ p = 0(l/e) (e — 0). 

In the proposition below, we call a net (s e ) e of complex numbers a slow-scale log- 
type net if there is a slow scale net (r e ) £ of real numbers, r e > 1, such that 

| SE |=0(log(r E )) (e-0). 

Proposition 3.4. In Theorem 13. II assume all log- type conditions to be replaced 
by slow-scale log-type estimates and, in addition, that (a £ ) E is of slow scale in each 
derivative. By the latter, we mean that for all j, k, I, we can find a slow scale net 
(r £ ) £ positive real numbers such that 

= 0(r £ ) (e - 0). 

Then / G 9^2 (X T ) and g £ 9^ 2 (R n ) implies u £ S£=,(Xr)- In particular, this is 
always true when the symbol of A is smooth (as in the Corollary above). 

Proof. Thanks to the explicit assumptions this is straightforward by an inspection 
of the proof of Theorem 13. II To be more precise, assume that we have a uniform 
e-growth, say £~ M , for the derivatives of / and g; i.e., for all k, a, we have 

\\d?d«f £ \\ L2 = 0(e- M ) as well as \\d^g £ \\ L2 = 0(e~ M ). 

Note that all constants involving a £ in the energy and Sobolev estimates throughout 
the proof yield only slow scale factors. (Observe again, that in the exponential 
factors in all energy estimates we only need a fixed finite order of derivatives, 
corresponding to k' n , l' n etc.) Thus, the same induction argument shows that we 
obtain for all k, a a certain slow scale net (r £ ) £ of positive real numbers, such that 

\\d?dSu £ \\ L2 =0(r £ e- M )=0(e- M - 1 ), 

whiches proves the assertion. □ 
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Remark 3.5. (i) The somewhat extensive slow-scale log-type conditions in the 
above proposition are by far not necessary for regularity, but are suited to make the 
energy estimates, with their exponential constants, directly applicable. We expect 
that these can be relaxed at least to plain slow scale conditions by appealing to 
pseudodifferential parametrix techniques (cf. (3ED- 

(ii) A slow-scale property of (a e ) e (in all derivatives) is implied, for example, by 
the log-type assumptions on (a £ ) £ if, in addition, only a S°°-type regularity of (a e ) e 
is assumed. This follows from ^] Proposition 1.6] and the fact that log(l/e) is a 
slow scale net. 



Appendix: Remainder term estimates 

We briefly outline a proof of l|lfi[) and verify the precise form of the estimate (|18f) ; it 
is an adaption of the reasoning in Ch.2, Sections 2-3]; thereby, we also recall the 
precise meaning of the oscillatory integral (|T7Jl . We may suppress the dependence 
of all symbols on the parameter t, since it will be clear that all steps in the process 
respect continuity (or smoothness) with respect to it and yield uniform bounds in 
all estimates when t varies in [0, T]. 

Let a{x, £) be a (smooth) symbol of order 1. The starting point is the following 
formula for the adjoint of a(x, D), e.g. valid for u G y(W l ) as iterated integral, 

a(x, D)*u{x) = J J e^ x ~ y ^ ri a(y,rj)u(y) dytTrj. 

Writing u(y) = J e lyi ctcj we obtain, now in the sense of oscillatory integrals, 

a(x,D)*u(x) = I e lx ^u(0 [[ e^ x - y ^a(y,^ + rj) dySr^ 
by which we define the symbol a* . Using Taylor expansion 



and (the oscillatory integral interpretation of) the Fourier identity 1 (9 r (a(., £,))){x) 
a(x, £) leads to 

a*(x,£) = a(x,£) + f ffv e l(x - y)ri d^a(y, £ + Or]) dyd-qdS. 



Noting that ■qe % ^ x ^ v ^' n — D y (e i ( x ~ v '' n ) and integrating by parts yields equations Ijl6(l 
and (fTTj) . We use the notation d x d^ — Y^j=i ^Xjd^j, re = J2j r j,o, and recall that 
(|17|l can be defined as the classical integral 



re(x,0 = // e-^>(l + \y\ 2 )- x (l-A 11 ) x (d x d i a(x + y^ + Qrj)) dydr,, 



where A > n/2, so that s e (x, £; y, rj) := (1 + |y| 2 )~ A (l - A,,) A (d x d(.a(x + y, £ + Wqj) 
is integrable. We have the estimate 

(33) \d^ S g(x^;y,rj)\< 

Cn, a ,0 ?i,+l + |a|,l+|/J|(°) (! + lyl)" 2A (l + K + 0V\ 2 r lal/2 , 
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where c n ,a,p is uniform in 6 6 [0, 1]. In order to prove l|18|) we have to estimate 

. =:Ji+/ 2 . 



M<l«l/2 kl>l£l/2 

For an upper bound of Zi we use and the implication \q\ < |£|/2 |£ + 077 1 > 
|£|/2 (when 6 e [0,1]) to find 

|il| < C n,a,/3,A92A+1+| Q |, 1+|/3| (°) ( X + ICI) - '" 1 , 

uniformly in 0. To estimate I2, we first use that e~ lyv — \r)\~ 2l (—A y ) l (e~' lvri ) and 
integrate by parts to obtain 

\h\<f J \vr 2l \(-\)'dpSsg(x^;y,v)\dyd v . 

r» |»,|>|£|/2 

We apply with (i replaced by 0+2lej (j = 1, . . . , n) to the integrand and arrive 
at 

|i'a|<c n)a , /J) « 1 A^x + i+|a|, 2 m+|/9|( a ) / l^r 2 '(l + le + ^| 2 )- |Q|/2 a-ry. 

hl>l?l/2 

By Peetre's inequality (1 + |£ + 6>7/| 2 )-H/2 < 2W/ 2 (1 + |£| 2 )-H/ 2 (1 + \0>q\ 2 )\ a U 2 , 
so that 

\h\ <<a,/J,U?2A+l+|a| j3 /+l+|/J|(o)(l + lel)" |a| / \v\~ ^ (1 + \H % ) W/3 % 

hl>l?l/2 

where the remaining integral is finite if 21 > n + \a\. 

Summing up, and combining the conditions 2A > n, 21 > n + \a\, we have shown 
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